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1. Introduction 

In these lectures I shall discuss the status of the determination of the three 
cosmological parameters which enter the Einstein equation and govern ge- 
ometry and evolution of space-time of the Universe: the Hubble constant 
Hq, the mass density parameter il. and the cosmological constant A. 

Among the three parameters, the Hubble constant is the dimensionfull 
quantity which sets the basic size and age of the Universe. The perennial 
effort to determine Hq dates back to Hubble (1925) and has a long history of 
disconcordance. Recent progress has done much to resolve the long-standing 
discrepancy concerning the extragalactic distance scale, but there are some 
newly revealed uncertainties in the distance scale within the Milky Way. 
The emphasis in this lecture is on discussion of these uncertainties. 

The mass density parameter directly determines the formation of cosmic 
structure. So, as our understanding of the cosmic structure formation is 
tightened, we should have a convergence of the Q parameter. An important 
test is to examine whether the parameter extracted from cosmic structure 
formation agrees with the value estimated in more direct ways. This gives 
an essential verification for the theory of structure formation. 

The third important parameter in the Friedmann universe is the cosmo- 
logical constant A. We now have some evidence for a non-zero A which, if 
confirmed, would have most profound implications for fundamental physics. 
This lecture will focus on the strength of this 'evidence'. 

We take the normalisation 



+ A = 1 



(1) 
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for the flat curvature, where A = A/3Hq with A the constant entering in 
the Einstein equation. The case with J7 = 1 and A = is referred to as the 
Einstein-de Sitter (EdS) universe. We often use distance modulus 

m - M = 5 log(ciL/10pc) (2) 

instead of the distance d^. For conciseness, we shah omit the units for the 
Hubble constant, (km s~^Mpc~^). 

After the Summer Institute there appeared several important papers on 
the distance scale. I try to incorporate these results in this article. 

2. The Hubble Constant 
2.1. HISTORICAL NOTE 

The global value of Hq has long been uncertain by a factor of two. Before 
1980 the dispute was basically between two schools: Sandage and collabora- 
tors had insisted on Hq = 50 (Sandage &; Tammann 1982); de Vaucouleurs 
and collaborators preferred a high value Hq = 90 — 100 (de Vaucouleurs 
1981). Conspicuous progress was brought by the discovery of an empirical 
but tight relationship between galaxy's luminosity and rotation velocity, 
known as the Tully-Fisher relation (Tully & Fisher 1977). The use of the 
Tully-Fisher relation has largely reduced subjective elements in the dis- 
tance work, and Hq = 80 — 90 has been derived from a straightforward 
reading of the Tully-Fisher relation. Representative of this work are the 
papers of Aaronson et al. (1986) and Pierce & Tully (1988). A doubt was 
whether the result was marred with the Malmquist bias — whether the 
sample selects preferentially bright galaxies, and hence the result was bi- 
ased towards a shorter distance (Kraan-Korteweg, Cameron Sz Tammann 
1988; Sandage 1993a). A related dispute was over the distance to the Virgo 
cluster, whether it is 16 Mpc or 22 Mpc: the different results depending on 
which sample one used. 

The next momentous advancement was seen in 1989—1990 when a few 
qualified distance indicators were discovered. One of them is a technique 
using planetary nebula luminosity function (PNLF), the shape of which 
looked universal (Jacoby et al. 1990a). Another important technique is the 
use of surface brightness fluctuations (SBF), utilizing the fact that the im- 
ages of distant galaxies show a smoother light distribution; while surface 
brightness does not depend on the distance, pixel-to-pixel fluctuations in a 
CCD camera decreases as d^^ (Tonry Sz Schneider 1988). They proposed 
that this smoothness can be a distance indicator if the stellar population 
is uniform. What was important is that the two completely independent 
methods predicted distances to individual galaxies in excellent agreement 
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with each other (CiarduUo, Jacoby & Tonry 1993). The PNLF/SBF dis- 
tance also agreed with the value from the Tully-Fisher relation, with a 
somewhat larger scatter. These new techniques, when calibrated with the 
distance to M31, yielded a value around Hq = 80 and the Virgo distance 
of 15 Mpc (For a review of the methods, see Jacoby et al. 1992). 

Around the same time the use of Type la supernovae (SNela) became 
popular (Tammann & Leibundgut 1990; Leibundgut & Pinto 1992; Branch 
Miller 1993). The principle is that the maximum brightness of SNIa is 
nearly constant, which can be used as an absolute standard candle. Arnett, 
Branch and Wheeler proposed that the maximum brightness is reliably 
calculable using models which are constrained from observations of released 
kinetic energy (Arnett, Branch & Wheeler 1985; Branch 1992). This led to 
Hq = 50 — 55, in agreement with the calibration based on the first Cepheid 
measurement of the nearest SNIa host galaxy using the pre-refurbished 
Hubble Space Telescope (HST) (Sandage et al. 1992). In the early nineties 
the discrepancy was dichotomous as whether Hq = 80 or 50. (see Fukugita, 
Hogan Sz Peebles 1993 for the status at that time; see also van den Bergh 
1989, 1994). 

The next major advancement was brought with the refurbishment mis- 
sion of HST, which enabled one to resolve Cepheids in galaxies as distant 
as 20 Mpc (1994). This secured the distance to the Virgo cluster and tight- 
ened the calibrations of the extragalactic distance indicators, resulting in 
Hq = (70 — 75) ± 10, 10% lower than the 'high value'. Another important 
contribution was the discovery that the maximum brightness of SNela varies 
from supernova to supernova, and that it correlates with the decline rate 
of brightness (Pskovskii 1984; Phillips 1993; Riess, Press &: Kirshner 1995; 
Hamuy et al. 1996a). This correction, combined with the direct calibration 
of the maximum brightness of several SNela with HST Cepheid observa- 
tions, raised the 'low value' of Hq to 65 ^^q, appreciably higher than 55. 
This seemed to resolve the long-standing controversy. 

All methods mentioned above use distance ladders and take the distance 
to Large Magellanic Clouds (LMC) to be 50 kpc {m — M = 18.5) as the zero 
point. Before 1997 few doubts were cast on the distance to LMC (TABLE 
1 shows a summary of the distance to LMC known as of 1997). With the 
exception of RR Lyr, the distance converged to m — M = 18.5 it 0.1, i.e., 
within 5% error, and the discrepency of the RR Lyr distance was blamed 
on its larger calibration error. It had been believed that the Hipparcos 
mission (ESA 1997) would secure the distance within MW and tighten the 
distance to LMC. To our surprise, the work using the Hipparcos catalogue 
revealed the contrary; the distance to LMC was more uncertain than we 
had thought, introducing new difficulties into the determination of Hq. In 
this connection, the age of the Universe turned out to be more uncertain 



4 



than it was believed. 

During the nineties, efforts have also been conducted to determine the 
Hubble constant without resorting to astronomical ladders. They are called 
'physical methods'. The advantage of the ladder is that the error of each 
ladder can be documented relatively easily, while the disadvantage is that 
these errors accumulate. Physical methods are free from the accumulation 
of errors, but on the other hand it is not easy to document the systematic 
errors. Therefore, the central problem is how to minimise the model depen- 
dence and document realistic systematic errors. Nearly ten years of effort 
has brought results that can be compared with the distances from ladders. 
The physical methods include the expansion photosphere model (EPM) for 
type II SNe (Schmidt, Kirshner & Eastman 1992) and gravitational lens- 
ing time delay (Refsdal 1964). Use of SNela maximum brightness was once 
taken to be a physical method (Branch 1992), but then 'degraded' to be a 
ladder, which however significantly enhanced its accuracy. 



TABLE 1. Distance to LMC as of 1997 



Method 


Ref 


Distance moduli 


Cepheid optical PL 


Feast & Walker 1987 


18.47±0.15 


Cepheid optical PL 


Madore & Freedman 1991 


18.50±0.10 


Cepheid IR PL 


Laney & Stobie 1994 


18.53±0.04 


Mira PL 


Feast & Walker 1987 


18.48±(0.06) 


SN1987A ring echo 


Panagia et al. 1991 


18.50±0.13 


SN1987A EPM 


Schmidt et al. 1992 


18.45±0.13 


RR Lyrae 


van den Bergh 1995 


18.23±0.04 



2.2. EXTRAGALACTIC DISTANCE SCALE 

The measurement of cosmological distances traditionally employs distance 
ladders (see Weinberg 1972). The most traditional ladders are shown in 
TABLE 2. The listings written in italic indicate new methods which cir- 
cumvent intermediate rungs. The most important milestone of the ladder is 
LMC at 50kpc [ra — M = 18.5). A distance indicator of particular historical 
importance (Hubble 1925) is the Cepheid period-luminosity (PL) relation, 
which is given a great confidence, but we note that it requires a few rungs 
of ladders to calibrate its zero point. 

Prior to the HST work there were only 4—5 galaxies with Cepheid dis- 
tances which could be used to calibrate secondary indicators. The reach of 
the ground-based Cepheid measurement is about 3 Mpc, which means that 
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TABLE 2. Traditional distance ladders 



Method 



Distance range typical targets 



Population I stars 

trigonometric or kinematic methods (ground) 
main sequence fitting (FG stars) Pop. I 
trigonometric method (Hipparcos) 
main sequence fitting (B stars) 
Cepheids [Population I] (ground) 
Cepheids [Population I\(HST) 
secondary (extragalactic) indicators 



<50 pc 

<200pc 

<500pc 

40pc-10kpc 

lkpc-3Mpc 

<30Mpc 

700kpc-100Mpc 



Hyades, nearby dwarfs 
Pleiades 

nearby open clusters 
open clusters 
LMC, M31, M81 
Virgo included 



Population II stars 

trigonometric method (Hipparcos) 
subdwarf main sequence fitting 
cluster RR Lyr 



<500pc 

lOOpc-lOkpc 

5kpc-100kpc 



nearby subdwarfs 

global clusters 

LMC, age determinations 



one cannot increase the number of calibrating galaxies from the ground. 
Pierce et al. (1994) could finally measure Cepheids in NGC 4571 in the 
Virgo cluster at 15 Mpc, but only with the best seeing conditions and dif- 
ficult observations. The refurbishment of HST achieved a sufficient power 
to resolve Cepheids at the Virgo cluster (Preedman et al. 1994). Now 28 
nearby spiral galaxies within 25 Mpc are given distances measured using 
the Cepheid PL relation (Ferrarese et al. 1999b). A typical random error is 
4-5% (0.08-0.10 mag), and the systematic error (from photometry) is 5% 
(0.1 mag) excluding the uncertainty of the LMC distance, to which the 
HST-Key Project(KP) group assigns 6.5% error (0.13 mag). The prime use 
of these galaxies is to calibrate secondary distance indicators which pene- 
trate into a sufficient depth that perturbations in the Hubble flow are small 
enough compared with the flow itself. 

Cepheids are Population I stars, so reside only in spiral galaxies. The 
calibration is therefore direct for TP and some SNela. Por early type galax- 
ies (fundamental plane or Z)„ — cr, and SEP) the calibration is not very 
tight; one must cither use some groups where both early and late galaxies 
coexist, or regard the bulges of spiral galaxies as belonging to the same 
class as early galaxies and avoid contaminations from discs. Additional ob- 
servations have been made for the galaxies that host SNcIa (Saha et al. 
1999). The results are summarised in TABLP 3. We include a few earlier 
SNIa results which employ a partial list of Cepheid calibrators. 

We accentuate the results with the two methods, SEP and SNela, in 
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TABLE 3. Hubble constant 



Secondary indicators 


Refs 


Hubble constant 


Tully-Fisher 


HST-KP (Sakai et al. 1999) 


71±4 ± 7 


Fundamental Plane 


HST-KP (Kelson ot al. 1999) 


78±8 ± 10 


SBF 


HST-KP (Ferrarese et al. 1999a) 


69±4 ± 6 


SBF 


Tonry et al. (1999) 


77±4 ± 7 


SNela 


Ricss ot al. (1995) 


67±7 


SNela 


Hamuy et al. (1996b) 


63±3 ± 3 


SNela 


Jha et al. (1999) 


64.4+^-^ 


SNela 


Suntzcff et al. (1999) 


65.6±1.8 


SNela 


HST-KP (Gibson et al. 1999) 


68±2 ± 5 


SNela 


Saha et al. (1999) 


60±2 


Summary (see text) 




(64 - 78) ± 7 



particular to those we underlined in the table. A cross correlation analysis 
showed that the relative distances agree well between SBF and others, 
including the Cepheid (Tonry et al. 1997; Freedman et al. 1997), and that 
it is probably the best secondary indicator presently available together with 
SNela; Also important is that there are now 300 galaxies measured with 
SBF, which are essential to make corrections for peculiar velocity flows for 
their < 4000 km s~^ sample (Tonry et al. 1999). (PNLF is an indicator 
of comparable quality, but it requires more expensive observations so that 
applications are rather limited; see Jacoby et al. 1996 for the recent work.) 
The final value of Tonry et al. from their / band survey is Hq = 77 zt 8, in 
which ±4 is allotted to uncertainties in the flow model and another ±4 to 
SBF calibration procedure in addition to the error of the Cepheid distance 
±6 (a quadrature sum is taken). There are a several other pieces of the 
SBF work to determine Hq, which generally result in Hq = 70 — 90 (e.g., 
Thomson ct al. 1997 using HST; Jensen et al. 1999 with K band; sec a 
review by Blakeslee et al. 1998). The new calibration made by the HST-KP 
group (Ferrarese et al. 1999a) would decrease Hq only by 2%. The difference 
in the flnal Hq between Tonry et al. and Ferrarese et al. comes from using 
different targets (the latter authors use only 4 clusters) and flow models. 

It is impressive that analyses of SNela Hubble diagram give virtually 
the same answer, even though the samples are all derived primarily from 
the Calan-Tololo sample of Hamuy et al. (1996b). A smaller Hq of Saha et 
al. (1997) basically reflects the absence of the the luminosity-decline rate 
correction, which pushes up Hq by 10%. The other notable difference is a 
slightly higher value of HST-KP (Gibson et al. 1999), who made a reanalysis 
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for all Cepheid observations performed by other groups and showed that 
their distances (to SN host galaxies) are all farther than would be derived 
from the HST-KP procedure. The average offset is as large as 0.16 mag 
(8%). This correction applies to all results other than HST-KP should we 
keep uniformity of the Cepheid data reduction. This is important especially 
when one compares the SN results with those from other secondary distance 
indicators, since the calibrations for the latter exclusively use HST-KP 
photometry. Taking the luminosity- decline rate correlation to be real and 
adopting Cepheid distance from the HST-KP data reduction, I adopt Hq = 
68 from SNela. 

We present two plots in Figure 1, (a) the estimates of maximum bright- 
ness of different authors and (b) the decline rate Ami5, the amount of the 
decrease in brightness over 15 days following maximum light, both as a 
function of metallicity [0/H]. The second plot shows how metallicity ef- 
fects are absorbed by the My — Ami5 relation and the first proves that 
there is little metallicity dependence in the corrected maximum brightness, 
though some scatter is seen among authors. 

Leaving out the uncertainty of the Cepheid distance, Hq from Tonry et 
al.'s SBF is 77±6, and that from SNela (HST-KP) is 68±4. The difference 
is 13%, and the two values overlap at Hq = 71. Allowing for individual 
two sigma errors, the overlap is in a range of Hq = 65 — 76. An additional 
uncertainty is 6% error {5Hq = ±4.5) from the Cepheid distance which is 
common to both, still excluding the uncertainty of the LMC distance. We 
may summarise Hq = 71 it 7 or 64—78 as our current standard, provided 
that LMC is at 50 kpc. All numbers in the table are within this range, 
except for the central value of Saha et al (1997). 

In passing, let us note that Hq = 75±15 (Preedman et al. 1997) obtained 
directly from the Cepheid galaxy sample agrees with the global value, im- 
plying that peculiar velocities are not so large even in a 10—20 Mpc region. 

This convergence is a great achievement, but keep in mind that the 
SNela results are still lower than those from other secondary indicators^ by 
10%. There are additional problems. First, all these analyses are based on 
a LMC distance modulus of m - M = 18.50 (Feast k Walker 1987; Madore 
&; Preedman 1991), which has recently been cast into doubt. In addition, 

remark is given to the TF distance. While Sakai et al. (HST-KP) derived Ho = 71± 
8 using GiovancUi et al.'s (1997) cluster sample, Tammann and collaborators (Tarmnann 
1999; Sandage & Tammann 1997) insist on a low value Hq = 53 — 56. Their cluster 
result (Federspiel et al. 1998) neglects the depth effect of the Virgo cluster: contrary to 
ellipticals, spiral galaxies are distributed elongating along the line of sight (Yasuda et al. 
1996). Hence identifying the centre of gravity of the spiral galaxy distribution with the 
true core leads to an offset. In fact the presence of substructure behind the Virgo core is 
confirmed with the Cepheid for NGC4649. Tammann et al.'s field result comes from the 
allocation of an unusually large dispersion to the TF relation, which largely amplifies the 
Malmquist bias. TuUy (1999) obtained Ho = 82± 16 (Tully et al. 1998). 
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Figure 1. (a) Maximum brightness of SNela (in the V band) adopted by different 
authors (see TABLE 3) as a function of [0/H] of host galaxies: solid circles, Gibson et al. 
(1999); solid triangles, Suntzeff et al. 1999; open circles, Jha et al. (1999); open triangles, 
Hamuy ct al. (1996b); open square, Saha et al. (1999). Note that Saha ct al.'s calibration 
is not necessarily brighter, which is mainly due to a different treatment of cxtiction. (b) 
The decline rate Amis measured in the B band (Phillips et al. 1999) as a function of 
[O/H] (Gibson et al. 1999). The slope of the curve is dAmis/diO/B.] ~ 0.28. 



metallicity effects could lead to systematic errors. Finally, should we derive 
the Hubble constant with the error of 10%, the problem of dust extinction 
could be an issue, and it is potentially coupled with the metallicity. We now 
consider these issues in greater detail. 



2.3. DISTANCE TO LMC 



The present status of the LMC distance is given in TABLE 4. The most 
traditional paths to the LMC distance follow the ladder shown in the upper 
half of TABLE 2. The Hipparcos satellite can measure a parallax down to 2 
milli arcsec (mas), corresponding to a distance of 500 pc (ESA 1997). It was 
a reasonable expectation that one could obtain the geometric distance to 
the Pleiades cluster, circumventing the main sequence fitting from nearby 
parallax stars to the Pleiades and thus securing the Galactic distance scale. 
Hipparcos observations have also opened a number of novel methods that 
can be used to estimate the distance to LMC. This and related activities, 
however, have actually brought confusions, rather than securing the dis- 
tance scale within the MW. We discuss several issues in order. 
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TABLE 4. Distance to LMC: Year 1997/1998 



MetVinH 


Ref 


T^lClfflTlPP TTioHnli 


Cepheid PL 


Feast & Catchpole 1997 


18.70 ± 0.10 




Paturel at al. 1997 


18.7 




Madore & Freedman 1998 


18.57 ± 0.11 




Luri et al. 1998 


18.29 ± 0.17 




Luri ct al. 1998 


18.21 ± 0.20 




(traditional) w/ new Pleiades 


18.26 


RR Lyrae (stat. para) 


Fernley et al. 1998 


18.31 ± 0.10 




Luri ct al. 1998 


18.37 ± 0.23 




Udalski 1998/Gould et al. 1998 


18.09 ± 0.16 


RR Lyrae (subdwarf) 


Reid 1997 


18.65 




Gratton et al. 1997 


18.60 ± 0.07 


Mira 


van Leeuwen et al. 1997 


18.54 ± 0.18 




Whitelock et al. 1997 


18.60 ± 0.18 


Red clump 


Udalski et al. 1998a 


18.08 ± 0.15 




btanek et al. 1998 


18.07 ± 0.04 




Cole 1998 


18.36 ± 0.17 


Eclipsing binaries 


Guinan et al. 1998 


18.30 ± 0.07 




(Udalski et al. 1998b) 


18.19 ± 0.13(?) 


SN 1987A Ring echo 


Gould & Uza 1998 


<18.37 ± 0.04 




Sonneborn et al. 1997 


18.43 ± 0.10 




Panagia et al. 1997 


18.58 ± 0.03 




Lundqvist & Sonnebornl997 


18.67 ± 0.08 


Cepheid PL 


Sekiguchi & Fukugita 1998 


18.10-18.60 




Sandage et al. 1999 


18.57±0.05 


Cepheid PL (BW method) 


Gieren et al. 1997 


18.49 ± 0.05 



2.3.1. "The Pleiades problem" 

The Pleiades cluster at 130 pc has been taken to be the first milestone of 
the distance work, since it has nearly solar abundance of heavy elements. 
This cluster is already too far to obtain a reliable parallax with the ground 
based observations, and its distance is estimated by tying it with nearby 
stars with solar metallicity employing main sequence fitting of FGK dwarfs 
(e.g., van Leeuwen 1983). The distance obtained this way agrees with an 
estimate via the Hyades, the nearest cluster to which geometric distance is 
available from the ground (Hanson 1980; van Altena et al. 1997), after a 
correction for large metallicity of the Hyades (VandenBerg &; Bridges 1984). 
It was then a natural exercise to confirm these estimates with a parallax 
measured by the Hipparcos. The result showed that the Pleiades distance is 



10 



shorter by 0.25 mag (12%) (van Leeuwen & Hansen-Ruiz 1997, Mermilliod 
et al 1997)! This is summarised in TABLE 5. 

Mermilhod et al.'s (1997) (see also dc Zeeuw et al. 1997) have shown 
that such a disagreement is seen not only for the Pleiades but also for 
other open clusters to some degree. A noteworthy example is that the lo- 
cus of the Praesepe ([Fe/H]=-)-0.095) agrees with that of the Coma Ber 
([Fe/H]=— 0.065) without mctallicity corrections, while we anticipate the 
former to be 0.25 mag brighter due to higher metellicity. 

This is a serious problem, since the disagreement means that either our 
understanding of FGK dwarfs, for which we have the best knowledge for 
stellar evolution, is incomplete, or the Hipparcos parallax contains system- 
atic errors (Pinsonneault et al. 1998; Narayanan & Gould 1999). The origin 
is not understood yet. 



TABLE 5. Pleiades distance summary 



Author/Method 




Distance modulus 


van Lccuwcn (1983) 
Linga (1987) 




5.57±0.08 
5.61 


Hyades (Ferryman et al. 1998) 
Pleiades— Hyades 
metallicity correction 


3.33±0.01 
2.52±0.05 
-0.22±0.03 


5.63±0.06 


van Leeuwen & Hansen-Ruiz 1997 
Mermilliod et al. 1997 
van Leeuwen 1999 




5.32±0.05 
5.33±0.06 
5.37±0.07 



2.3.2. Metallicity effects in the LMC Cepheid calibration 
The Cepheid distance to LMC is based on the calibration using open cluster 
Cepheids, the distances to which are estimated by B star main sequence 
fitting that ties to the Pleiades (Sandage & Tammann 1968, Caldwell, 1983, 
Feast & Walker 1987, Laney &; Stobie 1994). Metallicity has been measured 
for some of these calibrator Cepheids (Fry & Carney 1997). The residual 
of the PL fit shows a strong metallicity (Z) dependence. This means either 
the Cepheid PL relation suffers from a large Z effect, or the distances to 
open clusters contain significant Z-dependent errors (Sekiguchi & Fukugita 
1998). A correction for this effect changes the distance to LMC in either 
way, depending upon which interpretation is correct. 
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This metallicity dependence problem can be avoided if parallaxes are 
used to find the distances to calibrator Cepheids. Attempts were made 
(Feast & Catchpolc 1997; Luri et al. 1998; Madore & Freedman 1998) using 
field Cepheids in the Hipparcos catalogue. Unfortunately, Cepheid parallax 
data are so noisy (only 6 have errors less than 30%) that they do not allow 
calibrations tighter than ladders. Another skepticism is that 2/3 of Cepheids 
in the nearby sample (e.g., 14/26 in the Fcast-Catchpole sample) are known 
to have companion stars, which would disturb the parallax (Szabados 1997). 

2.3.3. Red clump 

The OGLE group revived the use of the red clump (He burning stage of Pop- 
ulation I stars) as a distance indicator. Paczynski &: Stanck (1998) showed 
that the / band luminosity of the red clumps depends little on metallicity 
(see. Cole 1998, however), and gave a calibration using the Hipparcos par- 
allax for nearby He burning stars. Udalski et al. (1998a) and Stanek et al. 
(1998) applied this to LMC, and obtained a distance modulus 18.1ib0.1, 
much shorter than those from other methods. This is a modern version of 
an analysis of Mateo & Hodge (1986), who reported 18.1±0.3. We should 
also recall that earlier analyses using MS fitting of OB stars resulted in 
a short distance of 18.2—18.3 (Schommer et al. 1984; Conti et al. 1986), 
though somewhat dismissed in the modern literature. 



2.3.4. Detached eclipsing binaries 

Detached double-spectroscopic eclipsing binaries provide us with a unique 
chance to obtain the distance in a scmi-gcomctric way out to LMC or even 
farther. From the information given by the light curve and velocity curve, 
one can solve for the orbital parameters and stellar radii (Andersen 1991, 
Paczynski 1997; Bell et al 1993 for an earlier application to LMC HV2226; 
Torres ct al. 1997 for an application to the Hyades). If surface brightness 
of the two stars is known from colour or spectrum, one can obtain the 
distance as d = {F/ f Y^^Ri where F and / are fluxes at the source and 
the observer and Ri is stellar radius. Guinan ct al. (1998) applied this 
method to HV2274 in LMC and derived m ~ M = 18.30 ± 0.07 with the 
aid of Kurucz' model atmosphere to estimate surface brightness from the 
spectrum. Udalski et al. (1998b) claimed that the extinction used is too 
small by an amount of AE{B—V) = 0.037 mag based on OGLE multicolour 
photometry. If we accept this correction the distance becomes 0.11 mag 
shorter, i.e., m — M = 18.19. 
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2.3.5. RR Lyr problems 

In the first approximation the luminosity of RR Lyr is constant, but in 
reahty it depends on metaUicity. The dependence is usually expressed as 

{Mv{RR Lyr)) = a[Fe/H] +b. (3) 

Much effort has been invested to determine a and b. The problem is again 
how to estimate the distance to RR Lyr. Unlike the case with Cepheids, 
there are no unique ladders for the calibration, and a variety of methods 
have been used, of which the best known is the Baade-Wesselink method. 
The calibration from the ground may be summarised as 

(My(RR Lyr)) = 0.2[Fe/H] + 1.04. (4) 

With this calibration we are led to the LMC distance of m — M ~ 18.3, as 
we saw in TABLE 1 above. 

The Hipparcos catalogue contains a number of field subdwarfs with 
parallax. This makes a ladder available to calibrate RR Lyr in globiilar 
clusters. Gratton et al. (1997) and Reid (1997) carried out this subdwarf 
fitting. Gratton et al. gave 

(My(RR Lyr)) = (0.22 ± 0.09) [Fe/H] + 0.76. (5) 

Their data are plotted in Figure 2, together with (4) and (5). Reid's result is 
also consistent. This zero point, being brighter by 0.3 mag (at [Fe/H]=~1.8) 
compared to (4), would bring the LMC distance to m — M = 18.5 — 18.6. 

There are a few analyses using the statistical parallax for field RR Lyr in 
the Hipparcos catalogue. Fernley et al. (1998) reported that their halo RR 
Lyr lie almost exactly on the curve of (4), rather than (5), and concluded 
a confirmation of the ground-based calibration. This is also endorsed by an 
analysis of Gould Sz Popowski (1998). 

The distance to eponymous RR Lyr was measured by Hipparcos. We 
see (Fig. 2) that RR Lyr itself does not fall on (5), but almost exactly on 
(4), although the error is fairly large. The uncertainties by 0.3 mag in the 
RR Lyr calibration translate to the LMC distance modulus 18.25—18.55. 

2.3.6. Conclusions on the LMC distance 

The distance to LMC is uncertain as much as 0.4 mag (20% in distance), 
ranging from 18.20 to 18.60. The results are rather bimodal around the 
two values close to the edges. A geometric method with SN1987A ring 
echo initiated by Panagia et al. 1991 does not differentiate between these 
two values: the data are too noisy and the result depends on the model 
of the light curve and emission lines that is adopted (Gould & Uza 1997; 
Sonneborn et al. 1997; see Fig. 8 of the latter literature for the data quality). 
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Figure 2. Calibrations of RR Lyr. The open points are taken from Gratton et al. (1997) 
with the dotted Une indicating (5). The soUd Une is the ground-based calibration (4). 
The solid point denotes the eponymous RR Lyr measured by the Hipparcos satellite. 

As we have seen in this section, recent observations with new techniques 
seem to tip the hst to the lower value. This is clearly a systematic effect, 
so that we cannot simply take an 'average of all observations'. Rather, we 
should leave both possibilities open. 

2.3.7. Age of the globular clusters 

The RR Lyr calibration is also crucial in the estimation of the age of glob- 
ular clusters, since the stellar age is proportional to the inverse of luminos- 
ity, i.e., inverse square of the distance. The modern evolution tracks of the 
main sequence agree reasonably well among authors. There are some dis- 
agreements in colours around the turn-off point, largely depending on the 
treatment of convection, but the luminosity is little affected (e.g., Renzini 
1991; Vandenberg et al. 1996, especially their Fig. 1). Absolute magnitude 
at the turn-off point My'-' of the main sequence is hence a good indicator 
of the age, as (Renzini 1991), 

logtg = -0.41 + 0.37Mf ° - QA3Y - 0.13[Fe/H] , (6) 

in units of Gyr, or 

logtg = -0.41 + (0.37a-0.13)[Fe/H]+0.37[(M^^-M^^)+6]]-0.43F, (7) 

if (3) is included. The difference of the magnitudes between the turn-off 
point and RR Lyr {Mf^ - M^^) varies little among clusters and is mea- 
sured to be 3.5±0.1 mag (Buonanno et al. 1989; see Chaboyer et al. 1996 
for a compilation) . The metallicity dependence of the cluster age disappears 
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if a = 0.35, i.e., the globular cluster formation is coeval (Sandage 1993b). 
Both (4) and (5), however, give a ~ 0.2,^ which indicates that metal-poor 
clusters appear older. 

The dichotomous calibrations of RR Lyr obviously affect the age of 
globular clusters. Another large uncertainty is whether the age-metallicity 
correlation is real, indicating metal-poor clusters formed earlier, or is merely 
due to a systematic error, with the formation of globular cluster being 
coeval. The possibilities are four-fold: 

h {m — M)i^y[C to (noncoeval) to(coeval) 
1.05 18.25 18Gyr 15Gyr 

0.75 18.55 14Gyr 12Gyr 

In addition there are ±10% errors from various sources (Renzini 1991; Bolte 
& Hogan 1995; VandenBerg et al. 1996; Chaboyer et al. 1996). Figure 3 

shows the age of various clusters from Gratton et al. (1997) and Chaboyer 
et al. (1998) both using the calibration close to (5). The [Fe/H] dependence 
is apparent. 

The claims of Gratton et al. (1997), Reid (1997) and Chaboyer et al. 
(1998) for young universe (11-12 Gyr) assume a coeval-formation inter- 
pretation together with the long RR Lyr calibration and take a mean of 
globular cluster ages. Three other possibilities, however, are not excluded. 

2.4. METALLICITY PROBLEMS WITH CEPHEIDS 

In most applications of the Cepheid PL relation, mctallicity effects are ne- 
glected, motivated by theoretical arguments that they will be very small. 
This results from double cancellations of the metallicity dependences be- 
tween core luminosity and atmosphere, as well as between the effects of the 
helium abundance and of heavier elements. The expected effect is (Stothers 
1988, Iben& Tuggle 1975; Chiosi et al. 1993) 

7;^ = dMx/d\Fe/B\ ~ ±0.05 dex mag~^ (8) 

for the A = y,/ pass bands. A new calculation of Sandage et al. (1999) 
gives |7a| < 0.1 for A = B,V,I. 

When one is concerned with a 10% systematic error in the cosmic dis- 
tance scale, the metallicity effect must be scrutinised. If it were as large 
as —0.5, say, the true Cepheid distance to normal spiral galaxies would be 

remark is made on a recent analysis of Kovacs & Jurcsik (1996), who obtained 
o < 0.19 from a model- independent approach using the Fourier coefficients of the light 
curves that correlate with the metal abundance. 
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Figure 3. Age of globular clusters as a function of [Fe/H]. The solid points are from 
Gratton et al. (1997) for three different stellar evolution models. The open points are 
from Chaboyer ct al. (1998). The soUd line shows (7) but offset by -0.06. The dashed 
line is 11.8 Gyr of Gratton ct al. 

longer by 10% relative to LMC ([O/H]=-0.4). The calibrator SNe used 
in earlier papers (SN1937C, SN1972E, 1981B and 1990N) all reside in low 
Z galaxies, but recent additions include SNe in high Z galaxies (1989B, 
and notably 1998bu), thus the sample spans a wider metallicity range 
(see Figure 1 above). There is now no relative difference in metallicity 
effects any more between the SBF and SNIa calibrator samples (the offset 
is A[0/H] < 0.1). Therefore, we cannot ascribe the difference in Hq to the 
metallicity effect of the Ccphcid PL relation: the effect slightly reduces Hq 
from both methods if the sign of 7 is negative. 

However, it is important to know the magnitude of 7. Observation- 
ally, Preedman & Madore (1990: FM) showed with the M31 data that the 
metallicity dependence is small {'Jbvri = —0.32 ± 0.21). Gould (1994), 
however, reanalysed the same data and concluded it to be as large as 
J = -0.88 ± 0.16. The EROS collaboration derived Jvi = -0.44 from 
a comparison between LMC and SMC (Beaulicu ct al. 1997). Kochanek 

(1997) suggested 7^/7 = —0.14 it 0.14 from a global fit of galaxies with 
Cepheid observations. The metallicity dependence for Galactic Cepheids 
discussed in section 2.3.2 corresponds to jvjhk ~ —2. Kennicutt et al. 

(1998) pointed out that the metallicity gradient of M31 used by Freedman 
&; Madore is a factor of three too large and argued that the above values 
should be jbvri = -0.94 ± 0.78 (FM) and -2.1 ± 1.1 (Gould). 

Kennicutt et al. (1998) derived from HST observations of two fields in 
MlOl that jvi = —0.24 it 0.16, which is the value currently adopted in 
metallicity dependence analyses of the HST-KP group. If this is the true 
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value, the effect on the distance scale is of the order of 5±3% {Hq gets 
smaller). I would emphasize, however, that independent confirmations are 
necessary for this 7 value, since the MlOl analysis is based only on V and 
/ bands, and the effect of extinction might not be completely disentangled. 

2.5. CROSS-CHECK OF THE CEPHEID DISTANCES 

2.5.1. Tests with geometric methods 

NGC4258 (M106) is a Seyfert 2 galaxy with H2O mascr emission from 
clouds orbiting around a black hole of mass 4 x IO^Mq located at the 
centre. Precise VLBA measurements of Doppler velocities show that the 
motion of the clouds is very close to Keplerian and is perturbed very little 
(Miyoshi et al. 1995). A complete determination is made for the orbital 
parameters including centripetal acceleration and a bulk proper motion of 
the emission system. This yields a geometric distance to NGC4258 to be 
7.2 ± 0.3 Mpc (Herrnstein et al. 1999). 

Maoz et al. (1999) measured the distance to NGC4258 using the conven- 
tional Cepheid PL relation, and gave 8.1ib0.4 Mpc with (m— M)lmc = 18.5. 
This distance is 13% longer than that from the maser measurement. The 
short LMC distance would bring the Cepheid distance in a perfect agree- 
ment with the geometric distance. This is, however, only one example, and 
it can merely be a statistical effect: the deviation is only twice the error, so 
it may happen with a chance probability of 5%. 

2.5.2. Further checks for M31 

A number of distance estimates are available for the nearest giant spiral 
M31, and they are shown in TABLE 6 (the underlined numbers are the 
zero point). Stanek &: Garnavich (1998) applied the red clump method to 
M31, and obtained m - M = 24.47 ± 0.06, which agrees with the M31 
modulus 24.44 ± 0.10 from the Cepheid based on the [m — M)lmc = 18.5 
calibration, whereas the same method gives 18.1 for the LMC distance. 
Namely, M31— LMC largely disagrees between the two. This discrepancy 
might be ascribed to a metallicity effect of either Cepheids or red clumps, or 
to the systematic error of either indicator. Mochejska et al. (1999) ascribed 
it to the error of the Cepheid distance from a crowded stellar population. 
On the other hand, tip of giant branch (TRGB) gives (m — M)]v[3i-lmc in 
agreement with the value derived from the Cepheid. The difference from 
PNLF is also consistent. The value from RR Lyr, however, is consistent 
with that from red clumps (the numbers in the Table is derived using (4) 
for the zero point). The results are dichotomous again. 
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TABLE 6. Relative distance of M31 to LMC 



method 


M31 


LMC 


M31-LMC 


refs. 


Cepheid 


24.44±0.10 


18.5 


5.94±0.10 


Ferrarese et al. 1999b 


red clump 


24.47±0.06 


18.07±0.04 


6.40±0.07 


Stanek & Garnavich 1998 


TRGB 


24.41±(0.19) 


18.5 


5.91±(0.19) 


Ferrarese et al. 1999b 


RR Lyr (B) 


24.50±(0.15) 


18.30 


6.20± 


Pritchet & vd Bergh 1989 


PNLF 


24.44 


18.56±0.18 


5.82±0.18 


Jacoby et al. 1990b 



2.6. PHYSICAL METHODS 

2.6.1. Expansion photosphere model (EPM) for Type II supemovae 

This is a variant of the Baade-Wessehnk method. If a supernova is a black 
body emitter one can calculate source brightness from temperature; the 
distance can then be estimated by comparing source brightness with the 
observed flux. In SNell atmosphere the flux is diluted due to electron scat- 
tering opacity. If this grcyncss is calculated source brightness can be in- 
ferred. Schmidt, Kirshner & Eastman (1992) developed this approach and 
obtained the distances to SNell in agreement with those from the ladder. 
The point is that EPM gives absolute distance without zero point calibra- 
tions. The Hubble constant they obtained is 73=b9 (Schmidt et al. 1994). 

A possible source of systematic errors is in the estimation of the tem- 
perature from the spectrum or colour. The SNell physics also might not be 
uniform, as we see occasional large scatters in a cross-correlation analysis. 

2.6.2. Gravitational lensing time delay 

When quasar image is split into two or more by gravitational lensing, we 

expect the time delay among images, arising from different path lengths and 
gravitational potentials among image positions. The time delay between 
images A and B takes the form 



Dls 
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(9) 



where 6 is the angular difference between the source and image, A(j) is the 
difference in the potential and Dij is the angular diameter distances. The 
time delay is observable if the source is variable, and can be used to infer 
Hq (Refsdal 1964). Crucial in this argument is a proper modelling of the 
mass distribution of the deflector. The DD/D factor depends on only 
weakly; its A dependence is even weaker. 
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The first case where Hq is derived is with the 0957+561 lens system. The 
deflector is complicated by the fact that a giant elliptical galaxy is embed- 
ded into a cluster. Falco, Gorcnstcin & Shapiro (1991) noted an ambiguity 
associated with a galaxy mass — cluster mass separation, which does not 
change any observed lens properties but affects the derived Hubble con- 
stant. One way to resolve this degeneracy is to use the velocity dispersion 
of the central galaxy (Falco et al. 1991; Grogin & Narayan 1996). Kundic 
et al. (1997b), having resolved a long-standing uncertainty about the time 
delay, obtained i^o = 64 it 13 employing the Grogin-Narayan model. Tonry 
& Franx (1998) revised it to 71+7 with their new velocity dispersion mea- 
surement near the central galaxy. More recently, Bernstein & Fischer (1999) 
searched a wider variety of models, also using weak lensing information to 
constrain the mass surface density of the cluster component, and concluded 
Hq = 77^24' th*^ large error representing uncertainties associated with the 
choice of models. 

The second example, PG1115+080, is again an unfortunate case. The 
deflector is elliptical galaxy embedded in a Hickson-type compact group of 
galaxies (Kundic et al. 1997a). Keeton Sz Kochanek (1997) and Courbin et 
al. (1997) derived (51 — 53) ± 15 from the time delay measured by Schcchtcr 
et al. (1997). Impey et al. (1998) examined the dependence of the derived 
Hq on the assumption for the dark matter distribution, and found it to vary 
from 44ib4 (corresponding to M/ L linearly increasing with the distance) to 
65 lb 5 (when M/L is constant over a large scale). The latter situation may 
sound strange, but it seems not too unusual for elliptical galaxies, a typical 
example being seen in NGC5128 (Peng et al. 1998). 

Recently, time delays have been measured for three more lenses, B0218+ 
357, B1608-b656 and PKS1830-211. B0218-F357 is a rather clean, isolated 
spiral galaxy lens, and Biggs (1999) derived Hq = 69 (the central value 
will be 74 if = 0.3) with a simple galaxy model of a singular isothermal 
ellipsoid. For B1608-I-656, Koopmans & Fassnacht (1999) obtained 64±7 
for n = 0.3 (59±7 for EdS). For PKS1830-211, they gave 75+}^ for EdS 
and 85^'^ for = 0.3 from the time delay measured by Lovell ct al. 
(1998). More work is clearly needed to exhaust the class of models, but these 
three lens systems seem considerably simpler than the first two examples. 
Koopmans Sz Fassnacht concluded 74+8 for low density cosmologies (69 ±7 
for EdS) from four (excluding the second) lensing systems using the simplest 
model of deflectors. It is encouraging to flnd a good agreement with the 
values from the ladder argument, though the current results from lenses 
are still less accurate than the ladder value. It would be important to ask 
whether Hq < 60 or > 80 is possible within a reasonable class of deflector 
models. 
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2.6.3. Zeldovich-Sunyaev effect 

The observation of the Zeldovich-Sunyaev (ZS) effect for clusters tells us 
about the cluster depth (times electron density), which, when combined 
with angular diameter (times electron density square) from X ray observa- 
tions, gives us the distance to the cluster provided that cluster is spherical 
(Birkinshaw et al. 1991, Myers et al. 1997). This is often taken as a physical 
method to measure Hq. I give little weight to this method in these lectures, 
since it is difficult to estimate the systematic errors. The currently avail- 
able results wildly vary from a cluster to a cluster. The most importnat is 
a bias towards elongation. None of the known clusters are quite spherical, 
and selection effects bias towards clusters elongated along the line of sight 
because of higher surface brightness. This may happen even if one uses a 
large sample. Additional systematics arise from the sensitivity of the ZS 
effect to the cluster envelopes; one must resort to a model to correct for 
this effect. 

2.6.4. Physical methods: summary 

Physical methods now yield the Hubble constant which can be compared 
with that from ladders. TABLE 7 presents a summary of Hq from the 
physical methods. However, effort is still needed to determine systematic 
errors associated with the use of specific methods. 



TABLE 7. Hubble constant from ladders and physical methods 

mctliod //() r(>f('r('iic(> 

ladders n±7(x0.95-1.15) 

physical: EPM 73±9 Schmidt et al. 1994 

physical: lensing (low O) 74±8 Koopmans & Fassnacht 1999 

(EdS) 69±7 

physical: ZS (54±14) Myers et al 1997 



2.7. CONCLUSIONS ON Hq 

The progress in determining the cxtragalactic distance scale has been dra- 
matic. The ladders yield values convergent within 10%, which is compared 
to a factor 1.6 disagreement in the early nineties. A new uncertainty, how- 
ever, becomes manifest in the Galactic distance scale: there is a 15—20% 
uncertainty in the distance to LMC. Therefore, we may summarise 



(10) 
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as a currently acceptable value of the Hubble constant. This agrees with 
that from HST-KP (Mould et al. 1999) up to the uncertainty from the 
LMC distance, though we followed a different path of argument. This allows 
Hq = 90 at the high end (if Tonry et al's SBF is weighted) and 60 at the low 
end (if the SNela results are weighted) . Note that Hq from both EPM and 
gravitational lensing arc consistent with the ladder value for (m— M)lmc = 
18.5. With the shorter LMC distance the overlap is marginal. 

The short LMC distance will also cause trouble for the i^o— age consis- 
tency. The LMC distance modulus of m — M = 18.25 would raise the lower 
limit of Hq to 72, and increase the lower limit of age from f«11.5 Gyr to 
~14.5 Gyr at the same time. There is then no solution for a A = universe. 
With a non-zero A, a unique solution is Hq ~ 72, $1 ~ 0.25, A ~ 0.75 with 
coeval globular cluster formation (see Figure 6 below). 

In the future it is likely that more effort will be expended for geometric 
methods. The great advantage is that it is free from errors arising from 
the chemical composition. In the surface brightness method, the chemical 
composition may still enter into the game, but its effect is tolerable and can 
even be reduced to a negligible level by using near infrared observations. 

Ultimately, gravitational wave observations could provide us with a 
novel method. For instance, for coalescing binary neutron stars the distance 
can be calculated as d ~ u^'^e^^T^^ , where e is metric perturbations, v is 
the frequency and r = i^/j> is a characteristic time of the collapse (Schutz 
1986). The position of the object may be difficult to infer, but there might 
be a gamma ray burst associated with the coalescence. 

3. The density parameter 

3.1. MODEL-INDEPENDENT DETERMINATIONS 

3.1.1. Luminosity density x [M / L) 

The mass density can be obtained by multiplying the luminosity den- 
sity with galaxy's average mass to light ratio (M/L). The local luminos- 
ity density, evaluated by integrating the luminosity function, is reason- 
ably well converged to Cb = (2.0 ± 0.4) x 10^ /iL© Mpc~^ from many 
observations. The M/Lb of galaxies generally increases with the scale. 
When the mass is integrated to !=a 100 kpc, a typical M/Lb is about 
(100 — 200)/i in solar units, and it may still increase outward (e.g., Faber 
& Gallagher 1979; Little & Tremaine 1987; Kochanek 1996; Bahcah et al. 
1995; Zaritsky et al. 1997). The virial radius in a spherical collapse model 
is r = 0.13 Mpc 9,-^-'^^[M /W^Mq^^I^^^^^^. If the dark matter distribution 
is isothermal within the virial radius, the value of M/Lb inside the virial 
radius is (150 — 400)/i for L* galaxies. This is about the value of M/Lb 
for groups and clusters, (150 — 500)/i. Multiplying the two values we get 
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O = 0.20 X 2^^. See also Fukugita, Hogan & Peebles (1998) for variants of 
this argument. 

Carlberg ct al. (1996, 1997a) tried to make the argument more quantita- 
tive using their cluster sample and a built-in field galaxy sample. They es- 
timated M/Lr ^ (210±60)/i for field galaxies from the cluster value (289 it 
50)/t. Their luminosity density of field galaxies is jO-r = (1.7±0.2) x 10^ HLq 
Mpc"3, and therefore = 0.19 ± 0.06. Note that M/Lb 1.4 x M/Lr in 
solar units for the respective pass bands. 

The important assumption for these calculations is the absence of co- 
pious matter outside the clusters. This is a question difficult to answer, 
but the observation of weak lensing around the clusters indicate that the 
distributions of dark mass and galaxies are similar at least in the vicinity 
of clusters (Tyson Sz Fischer 1995; Squires et al. 1996). 

Some attempts have also been made to estimate the mass on a superclus- 
ter scale. Small et al. (1998) inferred M/Lb — 560/i for the Corona Borearis 
supercluster, by applying the virial theorem (inspired by an N body sim- 
ulation). On the other hand. Kaiser et al. (1998) estimated M/Lb — 250 
from a mesurcmcnt of the gravitational shear of weak lensing caused by a 
supercluster MS0302+17 the result is not well convergent, but it seems 
unlikely that CI is larger than 0.5. 

3.1.2. Hq versus cosmic age 

For Ho > 60, the age is 10.9 Gyr for the EdS universe. This is too short. 
must be smaller than unity. If we take to > 11.5 Gyr n < 0.7. The limit is 
weak, but the significance is that EdS universe is nearly excluded. 

3.1.3. Type la supernova Hubble diagram 

The type la supernova Hubble diagram now reaches z ~ 0.4 — 0.8. It can 
be used to infer the mass density parameter and the cosmological constant. 
As we discuss later (section 4.1) the observation favours a low Q and a 
positive A. If we take their formal errors, < 0.1 is allowed only at three 
sigma for a zero A universe (Ricss et al. 1998; Perlmutter et al. 1999). A 
zero A open universe may not be excluded yet if some allowance is taken 
for systematic effects, but EdS geometry is far away from the observation. 
The best favoured value is approximately, 

J^«0.8A-0.4. (11) 

^They suggest f2 ~ 0.04 on the basis that only early-type galaxy population traces the 
mass distribution and the luminosity density is multiplied by the fraction of early-type 
galaxies (20%). It seems possible that late type galaxies reside in low density regions, 
causing only a small shear, which is buried in noise, and escaped from the measurement. 



22 



3.1.4. Baryon fractions 

A cluster is a virialised object with the coohng time scale longer than the 
dynamical time scale, and hence the physics is governed only by gravity (ex- 
cept for cooling flows in high density regions). The gas in clusters is shock 
heated to the virial temperature T ~ 7 x 10^(cr/1000km s~^)^ K, and thus 
emits X rays by thermal bremsstrahlung. From the luminosity and temper- 
ature of X rays one can infer the mass of the X ray emitting gas. It has 
been known that the gas amounts to a substantial fraction of the dynam- 
ical mass, which means that baryons reside more in the gas than in stars 
by an order of magnitude (Forman & Jones 1982). The argument was then 
elaborated by White et al. (1993b) based on ROSAT observations. From 19 
clusters White & Fabian (1995) obtained Mgas/Mgrav = 0.056^-^/^, where 
Mgrav is the dynamical mass. By requiring that the cluster baryon fraction 
agrees with in the field, we have = 0.066h^^^^r]iQ = 0.39(?7io/5), 

where rjio is the baryon to photon ratio in units of lO"^'^ and the last 
number assumes h = 0.7. 

An independent estimate is made from the Zeldovich-Sunyaev effect 
observed in clusters (Myers et al. 1997; Grego et al. 1999): Mgas/Mgrav = 
0.082/1-1 is translated to Q = 0.044/i-i??io = 0.31(??io/5). 

If we insert a probable value of the baryon to photon ratio from primor- 
dial nucleosynthesis calculations, r)io = 3 — 5, we have O = 0.2 — 0.4. 

3.1.5. Peculiar velocity - density relation 

This is one of the most traditional methods to estimate the cosmic mass 
density. The principles are spelled out by Peebles (1980). There are two 
basic tools depending on the scale. For small scales (r < 1 Mpc) the pertur- 
bations developed into a non-linear regime, and the statistical equilibrium 
argument is invoked for ensemble averages that the peculiar acceleration 
induced by a pair of galaxies is balanced by relative motions (cosmic virial 
theorem). For a large scale (r > 10 Mpc), where perturbations are still in 
a linear regime, the basic equation is 

V -v + Hon^-^S = (12) 

with 6 the density contrast. The contribution from a cosmological constant 
is negligible. The problem inherent in all arguments involving velocity is 
the uncertainty regarding the extent to which galaxies trace the mass dis- 
tribution (biasing), or how much mass is present far away from galaxies. 

Small-scale velocity fields: The status is summarised in Peebles (1999), 
where he has concluded f7(10kpc ^ r < IMpc) = 0.15 it 0.10 from the pair 
wise velocity dispersion (with samples excluding clusters) and the three 
point correlation function of galaxies via a statistical stability argument. 
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Bartlett k. Blanchard (1996) argued that it is possible to reconcile the ob- 
served velocity dispersion with ~ 1 if one assumes galactic halo extended 
beyond > SOOkpc. As Peebles (1999) argued, however, the halo is unlikely 
to be extended that much as indicated by the agreement of MW's mass at 
100-200kpc and the mass estimate for MW+M31 in the Local Group. 

Beyond a 10 Mpc scale, linear perturbation theory applies. An integral 
form of (12) for a spherical symmetric case {v/Hqv = ^^'^{6)/?)) applied to 
the Virgocentric flow gives O ~ 0.2 for v ~ 200 — 400km s~^ and {5) ~ 2, 
assuming no biasing (Davis and Peebles 1983). Recently, Tonry et al. (1999) 
argued that the peculiar velocity ascribed to Virgo cluster is only 140 km 
s~^, while the rest of the peculiar velocity flow is attributed to the Hyd-Ccn 
supercluster and the quadrupole field. For this case Q. ~ 0.06. We may have 
f2 ~ 1 only when half the mass is well outside the galaxies. 

Peebles (1995) argued that the configuration and kinematics of galaxies 
are grown following the least action principle from the nearly homogeneous 
primeval mass distribution. Applying this formalism to Local Group galax- 
ies, he inferred Q = 0.15 it 0.15. On the other hand, Branchini &: Carlberg 
(1994) and Dunn & Laflamme (1995) argued that this conclusion is not 
tenable if mass is distributed smoothly outside galaxies as in O = 1 CDM 
models. This seems, however, not very likely unless mass distribution is 
extended over 10 Mpc scale (Peebles 1999). 

Large scale velocity fields: There are a few methods to analyse the large- 
scale velocity fields based on (12). The direct use of (12) is a comparison 
of the density field derived from redshift surveys with measured peculiar 
velocities. Alternatively, one may use the density field reconstructed from 
observed velocity field for comparison with the actual density field, as in the 
POTENT programme (Dekel et al. 1990). A variant of the first method is to 
observe the anisotropy in redshift space (redshift distortion) (Kaiser 1987). 
As linear theory applies, Q. always appears in the combination [3 = Q^'^ /h 
where 6 is a linear biasing factor of galaxies against the mass distribution 
and can be inferred through non-linear effects. Much effort has been in- 
vested in such analyses (see e.g., Strauss & Willick 1995; Dekel et al. 1997; 
Hamilton 1998), but the results are still controversial. The value of 
derived from many analyses varies from 0.3 to 1.1, though we see a general 
trend to favour a high value. Notably, the most recent POTENT analysis using 
the Mark III compilation of velocities (Willick et al. 1997) indicates a high 
density universe O = 0.5 — 0.7, and O > 0.3 only at a 99% confidence level 
(Dekel et al. 1999). 

The difficulty is that one needs accurate information for velocity fields, 
for which an accurate estimate of the distances is crucial. Random errors 
of the distance indicators introduce large noise in the velocity field. This 
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seems particularly serious in the POTENT algorithm, in which the derivative 
V-v/fi^'^ and its square are numerically computed; this procedure enhances 
noise, especially for a small J7. The difficulty of inferring large scale velocity 
field may also be represented by the 'great attractor problem'. Lynden- 
Bell et al. (1988) found a large-scale velocity field towards the Hyd-Cen 
supercluster, but also argued that this supercluster is also moving towards 
the same direction attracted by a 'great (giant) attractor'. With Tonry 
et al.'s (1999) new estimate of the distance using SBF, this velocity field 
is modest, and Hyd-Cen itself serves as the great attractor that pulls the 
Virgo cluster, with the conclusion that CI is small. 

3.2. MODEL-DEPENDENT DETERMINATIONS 

The following derivations of the mass density parameter are based on the 
hierarchical clustering model of cosmic structure formation assuming the 
cold dark matter dominance. The extraction of Q is, therefore, indirect, 
but on the other hand, it is reasonable to appeal to such models since Q 
is the parameter that predominantly controls structure formation. Note 
that CDM model is the only model known today that successfully predicts 
widely different observations, yet there are no observations strong enough 
to refute its validity. We do not discuss results from cosmological models 
where physical processes other than gravity play a major role. 

3.2.1. Shape parameter of the transfer function 

The initial perturbations of the density fluctuation P{k) = \Sk\'^ ~ fc"' re- 
ceive a modification as P{k) = ~ k^T{k) as they grow, where T{k) 
is called the transfer function. Fluctuations of a small scale that enter the 
horizon in the radiation dominant epoch do not grow for a while, till the uni- 
verse becomes matter dominated. The transfer function T(k) thus damps 
for small scales as ~ k~^, whereas it stays close to unity for long- wave 
lengths. The transition region is controlled by a parameter k ~ 27r/ctcq, 
cteq being the horizon size at the time of matter-radiation equality, i.e., 
a characteristic length of 6.5{Qh)~^h~^ Mpc. The parameter T = de- 
termines the behaviour of the transfer function and is called the shape 
parameter. To give a sufficient power to several tens of Mpc, F must be as 
small as 0.2 (Efstathiou et al. 1990). This small value (F = 0.15 - 0.25) is 
supported by later analyses (e.g.. Peacock & Dodds 1994; Eke et al. 1998). 

3.2.2. Evolution of the rich cluster abundance 

The cluster abundance at z ^ requires the rms mass fluctuation as = 
{{6M/M))^/'^\^^sh-mpc to satisfy (White et al. 1993a; Eke et al. 1996; Pen 
1998; Viana & Liddle 1999; see also Henry & Arnaud 1991) 
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0-8 ~ Om-^-^ . (13) 

The evolution of the cluster abundence is sensitive to ag in early epochs 
of growth for a given mass; it is z ^ 0.3 for rich clusters. The rich cluster 
abundance at z ^ 0.3 — 1, when compared with that at a low z, determines 
both <T8 and fl (Oukbir & Blanchard 1992). Carlberg et al. (1997b) derived 
n = 0.4 ±0.2, and Bahcall & Fan (1998) obtained n = 0.2+°;^ correspond- 
ing to a slow growth of the abundance. On the other hand, Blanchard & 
Bartlett (1998) obtained ~ 1 from a more rapid growth. A high value 
is also claimed by Reichart et al. (1999), while Eke et al. (1998) reported 
CI = 0.43 zb 0.25 for an open, and ^ = 0.36 ± 0.25 for a flat universe. 

The controversy among authors arises from different estimates of the 
cluster mass at high z. This is a subtle effect, since the mass varies little 
over the range of relevant redshift, while the cluster number density evo- 
lution is sufficiently rapid at fixed mass (Pen 1998). At low z we have an 
established mass temperature relation, and the cluster mass is securely esti- 
mated (Henry & Arnaud 1991). At high z, however, such direct information 
is not available. Blanchard &; Bartlett and Eke et al. used mass tempera- 
ture relations as a function of z derived from hydrodynamic simulations. 
Reichart et al. used an extrapolated mass X-ray luminosity relation. Bah- 
call and Fan used more direct estimates of the cluster mass at higher z 
for three clusters. A change of a factor of two in the mass estimate would 
modify the conclusion. 

3.2.3. Cluster abundance versus the COBE normalisation 
There are a number of ways to infer cxg from galaxy clustering and peculiar 
velocity fields. The problem with the information from galaxy clustering 
is that it involves an unknown biasing factor, which hinders us from de- 
termining an accurate erg. The velocity data are susceptible to noise from 
the distance indicators. Therefore, the cluster abundance discussed above 
seems to give us a unique method to derive an accurate estimate of ag for a 
low z universe. Another place we can extract an accurate erg is the fluctu- 
ation power imprinted on cosmic microwave background radiation (CBR) 
anisotropies. Currently only the COBE observation (Bennett et al. 1996) 
gives sufficiently accurate cjg = (Tg{HQ,Vl,X,VLB, ■■■)■ Assuming the model 
transformation function, the matching of COBE erg with that from the 
cluster abundance gives a significant constraint on cosmological parame- 
ters 9, = n{Ho,X) (Efstathiou et al. 1992; Eke et al. 1996)Figure 4 shows 
allowed regions for two cases, open and fiat universes, assuming a flat per- 
turbation spectrum n = 1 and ignoring possible tensor perturbations. 

The transfer function is modified if n 7^ 1. The possible presence of the 
tensor perturbations in CBR anisotropies causes another uncertainty. The 
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Figure 4- Parameter regions allowed by matching the rms fluctuations from COBE with 
those from the cluster abundance. A flat spectrum (n = 1) is assumed and the tensor 
perturbations arc neglected. The lower band is for a flat universe, and the upper one for 
a universe with A = 0. 

COBE data alone say n being between 0.9 and 1.5 (Bennett et al. 1996), but 
the allowed range is narrowed to n = 0.9 — 1.2 if supplemented by smaller 
angular-scale CBR anisotropy data (Hancock et al. 1998; Lineweaver 1998; 

Efstathiou ct al. 1998; Tegmark 1999). The presence of the tensor mode 
would make the range of n more uncertain as well as it reduces the value 
of as- The limit of n when the tensor mode is maximally allowed is about 
< 1.3^. Notwithstanding these uncertainties, 17 > 0.5 is difficult to reconcile 
with the matching condition. On the other hand, a too small fl (^ 0.15) is 
not consistent with the cluster abundance. 

3.2.4. Power spectrum in nonlinear galaxy clustering 

Peacock (1997) argued that the power spectrum in a small scale region 
{k~^ < 3h~^ Mpc), where nonlinear effect is dominant, shows more power 
than is expected in Q = 1 cosmological models. He showed that the excess 
power is understood if the mass density is 0.3. 

3.2.5. CBR anisotropy harmonics 

The I distribution of the CBR harmonics depends on many cosmological 
parameters. Precise measurements of the harmonics will allow an accurate 
determination of the cosmological parameters up to geometrical degeneracy 
(Zaldarriaga et al. 1997; Efstathiou & Bond 1999; Eisenstein et al. 1999). At 

''in Tegmark's analysis n < 1.5 is quoted as an upper bound, but this is obtained 
by making Qb (and Ho) a free parameter. If one would fix the baryon abundance, the 
allowed range is narrower, n < 1.3. 
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present the data do not give any constraint on but on some combination 
of and A; so we defer the discussion to the next section. 

4. Cosmological constant 

Currently three tests yield useful results on the problem as to the existence 
of the cosmological constant: (i) the Hubble diagram for distant type la 
supernovae; (ii) gravitational lensing frequencies for quasars; (iii) position 
of the acoustic peak in the harmonics of CBR anisotropics. 

4.1. TYPE lA SUPERNOVA HUBBLE DIAGRAM 

The luminosity distance receives a cosmology dependent correction as z 
increases; in a way O pulls down di and A pushes it up. (In the first order 
of z the correction enters in the combination of go = — X, so this 
is often referred to as a qq test.) The discovery of two groups that distant 
supernovae are fainter than are expected from the local sample, even fainter 
than are expected for qq = 0, points to the presence of A > (Riess et al. 
1998; Schmidt et al. 1998; Perlmutter et al. 1999). 

The general difficulty with such a Hubble diagram analysis is that one 
has to differentiate among a few interesting cosmologies with small differ- 
ences of brightness. For instance, at z = 0.4 where many supernovae are 
observed, the difference is Am = 0.12 mag between (0,A) = (0.3,0.7) and 
(0,0), and Am = 0.22 from (0,0) to (1.0,0). Therefore, an accuracy of 
(^5%) must be attained including systematics to conclude the presence of 
A. On the other hand, there are a number of potential sources of errors: 

(i) K corrections evaluated by integrating spectrophotometric data that are 
dominated by many strong features; 

(ii) relative fluxes at the zero point (zero mag) across the colour bands; 

(iii) dust obscuration in a host galaxy; 

(iv) subtraction of light from host galaxies; 

(v) identification of the maximum brightness epoch, and estimates of the 
maximum brightness including a Ami5 correction; 

(vi) selection effects (for high z SNe); 

(vii) evolution effects. 

Except for (vii), for which wc cannot guess much^, the most important 
seems to be combined effects of (i), (ii) and (iii). It is not easy a task 
to reproduce a broad band flux by integrating over spectrophotometric 
data convoluted with filter response functions, especially when spectrum 

^Riess et al. (1999) showed that the rise time is different between low z and high z 
samples, indicating some evolution of the SNIa population. The effect on the cosmological 
parameter is not clear. 
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TABLE 8. Estimates of maximum brightness on SNe: 
1997 vs. 1999 from Perlmutter et al. (1997; 1999). 



SN 


1997 value 


1999 value 


difference 


SN1992bi 


(23.26±0.24) 


23.11±0.46 


(0.15) 


SN1994H 


22.08±0.11 


21.72±0.22 


0.36 


SN1994al 


22.79±0.27 


22.55±0.25 


0.24 


SN1994F 


(21.80±0.69) 


22.26±0.33 


(-0.58) 


SN1994am 


22.02±0.14 


22.26±0.20 


-0.24 


SN1994G 


22.36±0.35 


22.13±0.49 


0.23 


SN1994an 


22.01±0.33 


22.58±0.37 


-0.57 



Note: The numbers in the parentheses are not used in the final result of the 1997 paper. 



contains strong features. (Even for the spectrophotometric standard stars, 
the synthetic magnitude contains an error of 0.02—0.05 mag, especially 
when the colour band involves the Balmer or Paschen regions.) Whereas 
Perlmutter ct al. assigns 0.02 mag to the error of (i) [and (ii)], a comparison 
of the two values of estimated maximum brightness in their 1997 paper 
(Perlmutter et al. 1997, where they claimed evidence for a high fl universe) 
and the 1999 paper (TABLE 8) shows a general difficulty in the evaluation 
of the K correction (the difference dominantly comes from different K 
corrections). Schmidt et al. claim that their K correction errors are 0.03% 
mag. Dust obscuration (iii) is also an important source of errors, since the 
error of (i)+(ii) propagates to E{B — V) and is then amplified with the R 
factor. So a 0.02 mag error in colour results in a 0.06 mag error in Ay. 

We note that each SN datum contains ±0.2 mag (20%) error. The issue 
is whether this error is almost purely of random nature and systemtics are 
controlled to a level of <0.05. 



4.2. GRAVITATIONAL LENSING FREQUENCIES FOR QUASARS 
The gravitational lensing optical depth is given by 

'DolDls 



dT = flFH^{l+ZL 



Dos 



where F = {IGn^UgagHQ'^) , and is a magnification factor. The cosmo- 
logical factor in (14) is very sensitive to the cosmological constant, when 
it dominates (Fukugita & Turner 1991). F is the astrophysical factor that 
depends on the galaxy number density Ug and the mass distribution of 
galaxies, which is usually assumed to be a singular isothermal sphere with 
velocity dispersion ag. Figure 5 shows a typical calculation for the expected 
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Figure 5. Gravitational lensing frequencies as a function of A in a flat universe. The 
expecetd number is given for 504 quasars of the HST Snapshot Survey sample. The shade 
means the region within a ±50% uncertainty. The observed number is 5 (dashed hne). 

number of strong lenses for 504 quasars of the HST Snapshot Survey (Maoz 
et al. 1993) sample: the observed number is 5 (4 if 0957+561 is excluded). 
The curve shows a high sensitivity to A for A > 0.7, but in contrast a nearly 
fiat dependence for a lower A. It is likely that A > 0.8 is excluded. On the 
other hand, a more stringent limit is liable to be elusive. Fifty percent un- 
certainty in the F factor, say, would change largely a limit on, or a likely 
value of, A. 

In order to acquire information for a smaller A, an accurate estimate is 

essential for the F factor, which receives the following uncertainties in: (1) 
the luminosity density and the fraction of early-type galaxies (the lensing 
power of E and SO galaxies is much higher than that of spirals, and F is 
roughly proportional to the luminosity density of early- type galaxies); (2) 
(Tg-luminosity relation (Faber- Jackson relation); (3) the relation between 
C7(dark matter) and (7(star); (4) the model profile of dark haloes, specifically 
the validity of the singular isothermal sphere approximation (note that dark 
matter distributions seem more complicated in elliptical galaxies than in 
spiral galaxies, see Fukugita &: Peebles 1999); (5) the core radius which leads 
to a substantial reduction in dr; (6) selection effects of the observations; 
(7) dust obscuration; (8) evolution of early-type galaxies. 

There are continuous efforts for nearly a decade that have brought sub- 
stantial improvement in reducing these uncertainties (Maoz & Rix 1993; 
Kochanek 1996; Falco et al. 1998). Nevertheless, the issue (1) still remains 
as a cause of a large uncertainty. While the total luminosity density is 
known to an uncertainty of 20% or so, the fraction of early type galaxies 
is more uncertain. It varies from 0.20 to 0.41 depending on the literature. 
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Including other items, it is likely that an estimate of F has a 50% uncer- 
tainty. For the curve in Figure 5 a change of F by ±50% brings the most 
likely value of A to 0.75 or 0.2. 

Kochanek and collaborators have made detailed considerations on the 
above uncertainties, and carried out elaborate statistical analyses. In their 
latest publication they concluded A < 0.62 at 95% confidence level from 
an optical sample (Kochanek 1996). They took the fraction of early- type 
galaxies to be 0.44 and assigned a rather small la error. (The predicted 
frequency comes close to the upper envelope of Fig. 5, and the observed 
number of lenses in the HST sample is taken to be 4) . If one would adopt a 
smaller early-type fraction, the limit is immediately loosened by a substan- 
tial amount. Since the uncertainty is dominated by systematics rather than 
statistical, it seems dangerous to give significance to statistics. Statistical 
significance depends on artificial elements as to what are assumed in the 
input. A similar comment also applies to the recent work claiming for a 
positive A (Chiba k, Yoshii 1997; Cheng & Krauss 1998). I would conclude 
a conservative limit being A < 0.8. 

4.3. HARMONICS OF GBR ANISOTROPIES 

This is a topic discussed repeatedly in this Summer Institute, so I will only 
briefly mention it for completeness. The positions of the acoustic peaks are 
particularly sensitive to Q, and A, and even low accuracy data available at 
present lead to a meaningful constraint on a combination of Q, and A. 

The first acoustic peak appears at = -7r(thc distance to the last- 
scattering surface)/(the sound horizon) (Hu & Sugiyama 1995). Its position 
ii is approximated as 

for the parameter range that concerns us. This means that the position of 
the acoustic peak is about i ~ 220 if J7 -I- A = 1, but it shifts to a high ^ 
as if A = 0. On the other hand, there is little power to determine 

ri separately from A, unless full information of is used. The harmonics 
Q measured at small angles revealed the acoustic peak (Scott et al. 1996), 
and its position favours a universe not far from flat (Hancock ct al. 1998). 
More exhaustive analyses of Lineweaver (1998) and Efstathiou ct al. (1999) 
show a limit Q, + A/2 > 0.52 (Icr). (The contours of the confidence level fall 
approximately on the curve given by (15) with ^i=constant.) This means 
that a zero A universe is already marginal, when combined with Q, from 
other arguments. If a fiat universe is chosen from CBR, a non-zero A will 
be compelling. 
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TABLE 9. Summary of Q and A. 



method Qo A? model used? 

Ho vs to < 0.7 

luminosity density +M/L 0.1-0.4 

cluster baryon fraction 0.15-0.35 

SNela Hubble diagram < 0.3 A « 0.7 

small-scale velocity field (summary) 0.2 ± 0.15 

(pairwisc velocity) 0.15 ± 0.1 

(Local Group kinematics) 0.15 ± 0.15 

(Virgocentric flow) 0.2 ± 0.2 

large-scale vel field 0.2—1 

cluster evolution (low Q sol'n) 0.2 j yes 

(high n sol'n) '^l yes 

COBE-cluster matching 0.35-0.45 (if A = 0) yes 

0.20-0.40 (if A ^ 0) yes 

shape parameter F 0.2 — 0.4 yes 

GBR acoustic peak free (if flat) > 1 — 2f2 yes 

> 0.5 (if open) yes 
gravitational lensing A < 0.8 

summary 0.15—0.45 (if open) 

0.2-0.4 (if flat) 

0.6-0.7(?) 



5. Conclusions 

The status of and A is summarised in TABLE 9. We have a reasonable 
convergence of the parameter towards a low value O, = 0.15 — 0.4. The 
convergence of U is significantly better with the presence of the cosmological 
constant that makes the universe flat. Particularly encouraging is that the 
parameters derived with the aid of structure formation models agree with 
each other. This is taken to be an important test for the cosmological model, 
just as in particle physics when many different phenomena are reduced to 
a few convergent parameters to test the model. There are yet a still highly 
discrepant results on ii, but it is not too difficult to speculate their origins. 
On the other hand, the current 'low Q,' means the values that vary almost 
by a factor of three and effort is needed to make these converge. 

The cosmological constant has been an anathema over many years be- 
cause of our ignorance of any mechanism that could give rise to a very small 
vacuum energy of (3 meV)^, and neither can we understand a zero cosmo- 
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Figure 6. Consistent parameter ranges in the Hq — Q, space for (a) a flat universe and (b) 
an open universe. A is the range of the Hubble constant when (m — M)lmc = 18.5. B or C 
is allowed only when the LMC distance is shorter by 0.3 mag, or longer by 0.1 mag. Note 
in panel (a) that most of the range of B is forbidden by the compatibility of age and Ho 
that are simultaneously driven by the RR Lyr calibration (section 2.7). Also note that the 
age range between wll.S Gyr and Ril4 Gyr is possible only with the interpretation that 
globular cluster formation is coevel (section 2.4). The most naturally- looking parameter 
region is given a thick shade. 



logical constant. In mid-nineties the atmosphere was changing in favour for 
a non-zero A. The prime motivation was the Hubble constant— age prob- 
lem, but the introduction of a non-zero A was helpful in many respects. One 
theoretical motivation was to satisfy flatness which is expected in inflation- 
ary scenarios (Peebles 1984). Ostrikcr & Stcinhardt (1995) proclaimed a 
'cosmic concordance' with a flat universe mildly dominated by A. By 1997, 
only one observation contradicted with the presence of a moderate value 
of A; this was the SNela Hubble diagram presented by the Supernova Cos- 
mology Project (Pcrlmutter et al. 1997); see Fukugita 1997. In the next two 
years the situation changed. Two groups analysing SNela Hubble diagram, 
including the Supernova Cosmology Project, now claim a low ft and a pos- 
itive A. On the other hand, the Hubble constant— age problem became less 
severe due to our cognition of larger uncertainties, especially in the age 
estimate. The indications from SNela Hubble diagram are very interest- 
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ing and important, but the conclusions are susceptible to small systematic 
effects. They should be taken with caution. We should perhaps wait for 
small-scale CBR anisotropy observations to confirm a nearly fiat universe 
before concluding the presence of A. 

In these lectures we have not considered classical tests, number counts, 
angular-size redshift relations, and magnitude-redshift relations of galaxies 
(Sandage 1961; 1988), in those testing for J7 and A. Unlike clusters or 
large scale structure, where no physics other than gravity plays a role, the 
evolution of galaxies is compounded by rich physics. Unless we understand 
their astrophysics, these objects cannot be used as testing candles. It has 
been known that galaxy number counts is understood more naturally with 
a low matter density universe under the assumption that the number of 
galaxies are conserved, but it is possible to predict the correct counts with 
an O = 1 model where galaxies form through hierarchical merging, by 
tuning parameters that control physics (Cole et al. 1994; Kauffmann et 
al. 1994). It is important to work out whether the model works for any 
cosmological parameters or it works only for a restricted parameter range. 
This does not help much to extract the cosmological parameters, but it can 
falsify the model itself. 

We have seen impressive progress in the determination of the Hubble 
constant. The old discrepancy is basically solved. On the other hand, a 
new uncertainty emerged in more local distance scales. The most pressing 
issue is to settle the value of the distance to LMC. There are also a few 
issues to be worked out should one try to determine Hq to an accuracy of 
a 10% error or less. They include understanding of metallicity effects and 
interstellar extinction. The future effort will give more weight to geometric 
or semi-geometric methods. From the view point of observations the work 
will go to infra-red colour bands to minimise these problems. 

In conclusion, I present in Figure 6 allowed ranges of Hq and f2 (and A) 
for the case of (a) flat and (b) open universes. With the flat case we cut the 
lower limit of O at 0.2 due to a strong constraint from lensing. An ample 
amount of parameter space is allowed for a flat universe. A high value of 
Hq > 82, which would be driven only by a short LMC distance, is excluded 
by consistency with the age of globular clusters as noted earlier. Therefore, 
we are led to the range i^o — 60 — 82 from the consistency conditions. 
For an open universe the coeval-formation interpretation is compelling for 
globular clusters, or else no region is allowed. The allowed Hq is limited to 
60 — 70. No solution is available if LMC takes a short distance. 

I would like to thank Rob Crittenden for his careful reading and many 
useful suggestions on this manuscript. This work is supported in part by 
Grant in Aid of the Ministry of Education in Tokyo and Raymond and 
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